Novel superhard sp 3 -carbon allotropes with odd and even ring topologies 
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Four novel sp 3 -carbon allotropes with 6, 8 and 16 atoms per primitive cell have been derived using 
a combination of met adynamics simulations and topological scan. A novel chiral orthorhombic 
phase oC16 (C222i) was found to be harder than monoclinic M-carbon and shows remarkable 
stability in the high pressure range. A second orthorhombic phase of Cmmm symmetry, by ~0.028 
eV/atom energetically lower than W-Carbon, can be formed from graphite at ^9GPa. In general, 
the mechanical response under pressure was found to depend on the structure topology, which 
reflects the way rings are formed from an initial graphene layer stacking. 



I. INTRODUCTION 

The quest for novel carbon materials with improved 
mechanical properties and tailored optical gap is a topic 
of high priority. Engineering new properties is tightly 
connected with the ability to predict crystal structures, 
which remains a crucial issue in both basic solid state 
research and modern materials science 1 . 

In the effort of anticipating superior materials for catal- 
ysis, hydrogen storage and gas segregation, structure pre- 
diction stands out for its capacity to efficiently indicate 
viable technological target compounds. The challenge 
consists in identifying metastable modifications that can 
exhibit interesting physical and chemical properties. 

Compression of graphite at high pressure and tempera- 
tures produces diamond^. Graphite cold compression on 
the contrary produces a hard and transparent product, 
different either from cubic or hexagonal diamond^ - —, but 
not fully characterized so far. Many recent studies deal 
with the nature of this metastable product. Several en- 
ergetically competing carbon phases were proposed (W- 
and M-carbon^, bet C4) as plausible structure solutions, 
based on estimating transition pressures, goodness of fit 
of X-ray diffraction data and band gapsii - — . The intrin- 
sic problematic of stacking faults in the pristine graphite, 
and Raman evidence of amorphisation suggest a mixture 
of different phases in the compressed material. The two 
energetically most preferable candidates so far (M- and 
W-carbon) can be described (in terms of topology) as 
corrugated graphene sheets interconnected by an alter- 
nating sequence of odd rings (pentagons and heptagons) 
fused into a 5+7 pattern. This odd-ring topology for- 
mally results from connecting puckered graphene layers 
aligned in a particular way. On the other hand, further 
compressing M- or W-carbon can produce different dia- 
monds poly types. Therefore, a larger variety of interme- 
diate hard structures can in principle be expected. 

In this Communication, we further unfold the struc- 
tural diversity of sp 3 -carbon phases. We base our 
approach on metadynamics simulations of structural 
transformations^^ and topological enumeration to ef- 
ficiently scan configuration space. We report on ener- 
getic, mechanical and electronic properties of four new 
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FIG. 1. Crystal structures of novel carbon phases. 0CI6- 
I and mC32 are characterized by 5+7 odd-odd ring pattern. 
mC12 represents a distinct 5+8 odd-even ring topology, while 
oC 16-11 contains even rings only, 4+6+8. 



tetrahedral carbon phases, and insist on the different un- 
derlying graphitic pattern connected with the formation 
of a particular topology. We show how distinct topolo- 
gies with 5+7 (odd-odd), but also 5+8 (odd-even) and 
4+6+8 (even-even-even) ring pattern can do for different 
mechanical responses. 



II. METHODS 

Efficient theoretical approaches to hypothetical car- 
bon modifications, based e.g. on random techniques, 
genetic (evolutionary) algorithms, or accelerated molecu- 
lar dynamics^ result in important discoveries supporting 
experiments^^. In some approaches the use of graph 
theoretical methods 19 represents a means of increasing 
the sampling efficiency of carbon configuration. It was 
indeed a graph-theoretical approach that allowed to de- 
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rive all possible sp 3 -carbon allotropes with 4 atoms per 
cell (including the recently rediscovered bet C^)^. Meta- 
dynamics on the other hand explores the energy land- 
scape along collective reaction coordinates, which in case 
of high pressure polymorphs is represened by the simu- 
lation box itself. While metadynamics does not require 
prior knowledge of the energy landscape under investiga- 
tion, its sampling efficiency improves on combining many 
independent runs started from different initial configura- 
tions. Additionally, the number of atoms per simulation 
box is critical for capturing a particular atomic configura- 
tion. Diamond and lonsdaleite are important metastable 
forms of carbon. They can appear in the same metady- 
namics run only if the number of atoms in the box is at 
least 4 and multiple thereof. Similarly, including a mini- 
mum of 3 atoms (or multiple thereof) is sufficient to find 
a dense carbon with quartz topology, recently suggested 
from evolutionary algorithms 20 . 

To systematically include known and find novel car- 
bon forms, metadynamics runs were performed on sim- 
ulation boxes comprising three, four, six, eight, twelve 
and sixteen carbon atoms, respectively. A similar ap- 
proach has been shown to work well in connection with 
plain MD to search for ice phase o 21 i 22 . Quasi-random 
four-connected nets were used as starting configurations. 
We note here that the application of metadynamics in 
this case is slightly different from its typical use for sim- 
ulation of crystal-crystal structural phase transitions, as 
recently reviewed in Refs. 23,24 . While in both cases the 
simulation cell is used as order parameter, here one in- 
stead starts from a disordered configuration in a small 
cell and searches for low-energy configurations represent- 
ing crystalline structures with a given number of atoms 
in the unit cell. Each run was typically consisted of 
about 25000 metasteps. Within each metastep MD was 
performed in the NVT ensemble for at least 0.5 ps at 
300 K. In these preliminary scans the tight binding Ter- 
soff potential 25 was used, which ensured rapid and reli- 
able structure evolution thanks to its good description 
of sp 2 /sp 3 carbons. Molecular dynamics in the NVT- 
ensemble was performed with the CP2K cod o 26 i 27 . Struc- 
ture diversity was judged by calculating vertex symbols, 
which contain information on all shortest rings meet- 
ing at each atom, and coordination sequences, as imple- 
mented in the TOPOS package 28 . Both topological de- 
scriptors are widely used e.g. for the topological charac- 
terization of zeolites 28 . Metadynamics trajectories con- 
tained many foam-like structures with mixed sp 2 -sp 3 car- 
bon atoms and a few sp 3 -allotropes. In case of new tetra- 
hedral structure types, inferred from their topology, ideal 
space groups and asymmetric units were identified with 
the Gavrog Systre package 29 . In a subsequent set of runs, 
candidate structures were studied with respect to their 
transformability into diamond by metadynamics simula- 
tions using SIESTA 30 ' 31 as DFT/MD layer. 

In the initial metadynamics runs the choice of large 
pressure values is less critical. It is rather the number 
of atoms in the simulation box that decides whether a 



particular topology can be visited at all within a sin- 
gle metadynamics run. The structures presented in the 
following were harvested from metadynamics runs per- 
formed at and 5 GPa, with 6, 8 and 16 atoms in the 
simulations box. 

On idealized structures variable cell conjugate-gradient 
relaxation was performed within density functional the- 
ory (GGA, PBE) as implemented in the SIESTA 
package 30,31 . Electronic states were expanded by a 
double- basis set with polarization functions (DZP). Core 
states (Is 2 ) were described by norm-conserving Troul- 
lierMartins pseudopotentials 32 . The charge density was 
represented on a real-space grid with an energy cutoff 
of 200 Ry. Forces were relaxed to less than 0.01 eV/A. 
Convergence with respect to the number of k-points was 
carefully checked. 



III. RESULTS AND DISCUSSION 

Small boxes of 2 and 3 atoms expectedly produced cu- 
bic diamond and quartz, respectively. With four atoms 
both cubic and hexagonal diamond (Lonsdaleite) were 
collected. From 6, 8 and 16 atoms metadynamics three 
novel structures were found, two monoclinic (mC12 and 
mC32, Fig. [1]) and one orthorhombic (0CI6-I, Fig. [1]). 
From further propagating 0CI6-I in metadynamics runs 
at 100 GPa, oC 16-11 (Cmmm) was found. Their symme- 
tries and structural parameters are summarized in Ta- 
ble HI All phases correspond to a stacking of corrugated 
graphene layers interconnected by alternating sequence 
of pentagons and heptagons (0CI6-I and mC32, Fig. [1]), 
like for M- and W-carbon. Alternatively, pentagons and 
octagons, or square and octagons can also be placed be- 
tween puckered graphitic layers like it is realized in mC12 
and 0CI6-II, respectively (Fig. [I]). Physical properties of 
the new allotropes are compared with those of known 
structures in Table [TTJ In terms of volume per atom, 
0CI6-I is the densest, hardest structure, closely followed 
by 0CI6-II and mC12. With a calculated band gap of 
4.5 eV, 0CI6-I is also the structure closest to diamond. 

The stability of different carbon phases in a wide pres- 
sure range is presented in Fig. [2j At elevated pressures, 
the novel allotropes become more stable than graphite, 
Fig-Gb- The transition pressures are similar for the mC32 
and 0CI6-I structures (19.7 and 23.4 GPa, respectively) 
and much higher (by 10 GPa) for mC12. mC12 and 
particularly 0CI6-I are stabilized upon increasing pres- 
sure. Furthermore, the stability of 0CI6-I remains ba- 
sically constant (up to 400 GPa) whereas M- and W- 
carbon rapidly become energetically unfavourable above 
100 GPa. 

Bulk moduli (Bo) were obtained by fitting total en- 
ergy as a function of volume to the third-order Birch- 
Murnaghan equation of state (Table [II]). Strickingly, 
0CI6-I is harder than M- and W-carbon, although less 
stable below 129.2 GPa (Fig. EJ. On the contrary, 0CI6- 
II features a lower enthalpy, its gap is nonetheless inter- 
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TABLE I. Crystal structure information for novel carbon phases at GPa 



Pearson 


Space Group 


Wyckoff 




y 




Symbol 


Cell 


position 


X 


z 


mC12 


C2/c 

a — 0.4Z4Z, D — o.OZlo, 


4e 





0.80280 


3/4 




c= 3.7012; /3=138.96° 


8f 


0.84662 


0.91988 


0.95940 


mC32 


C2/m 


8j 


0.46444 


0.68220 


0.12680 




a=9.7242; b=4.2932; 


8j 


0.94998 


0.68122 


0.59997 




c=4.8617; /3=103.96° 


8j 


0.30907 


0.68472 


0.43555 






8j 


0.18908 


0.68609 


0.87688 


0CI6-I 


C222i 


4a 


0.43209 


1/2 







a=6.6698; b=5.5609; c=2.5119 


4b 


1/2 


0.08196 


1/4 






8c 


0.81701 


0.76297 


0.11960 


0CI6-II 


Cmmm 


8p 


0.66672 


0.68505 


0.0 




a=8.8134; b=4.2743; c=2.5281 


8q 


0.58903 


0.81586 


1/2 



8.9 GPa 19.7 GPa 23.4 GPa 27.3 GPa 
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FIG. 2. (a) Enthalpies (relative to graphite) of different car- 
bon allotropes; (b) enthalpies of certain sp 3 -allotropes (rela- 
tive to diamond) in the high pressure range. The colours are 
the same as in (a). 



mediate between lonsdaleite and bet C4, which is struc- 
turally also the case. By inspection of Fig. [TJ motifs of 
lonsdaleite and bet C4 can easily be recognized. 

The electronic band structures of mC12, 0CI6 (I and 
II) and mC32 at 50 GPa are shown in Fig. [3] The struc- 
tures are insulating with indirect band gaps in the range 
2.8-4.5 eV. The gaps do not depend on the pressure up 
to 50 GPa. All the gaps are smaller than in diamond, 
but similar to those of M- and bet-carbon. 

Phonon dispersion curves were calculated within a 
pressure range up to 100 GPa. No imaginary frequen- 
cies were observed throughout the whole Brillouin zone, 
confirming the dynamical stability of the intermediate 
sp 3 structures (Fig.[3j). Isothermic-isobaric molecular dy- 
namics simulations (300 K, 1 atm, 3 ps) also confirmed 
the stability of the new phases. 

Fig. [4] shows the relation between the discovered struc- 
tures and the graphene layer stackings they are deriverd 
from. This information can be obtained by deconstruct- 
ing the structures and looking for graphitic layers within 
the lattices. For the novel structures, the matching we 
are presenting is supported by the met adynamics runs, 
where a (fully or partially) graphitic structure is a typ- 
icall precursor of the sp 3 phases, along the simulation 
time coordinate. 

In general, we notice that hardness and band gap are 
diver sly distributed among the phases. In the effort of 
providing an answer to the outstanding question of hard 
and transparent sp 3 carbon, oC 16-11 and 0CI6-I appear 
as better candidates as hitherto suggested, the former 
particularly for its stability, for a really transparent band 
gap and hardness the latter. In Fig. [5] we present the 
simulated XRD patterns of 0CI6-I and oC 16-11. With 
reference to the experimental pattern^, the relevant re- 
gions between 8.5-10° as well as 14.5-17° are similarly 
populated. Intermediate peaks between 10-14° can bet- 
ter distinguish between the two structures, are however 
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TABLE II. Calculated equilibrium volume (Vo), bulk modulus (Bo) and band gaps (E g ). All values refer to zero pressure. 



Structure 


iviexnoa 


VO (A J 


DO ^O-Taj 




ri(Lrpaj_ 




This work 


^ 7Q 

0. ( y 


/1 9/1 9 

4Z4.Z 


/I 1 Q 

4. iy 


87 


Diamond 


PBE 20 


5.70 


431.1 


4.2 




mC12 


This work 


5.91 


399.5 


2.82 


84.4 


n pi a t 
Ov^lD-1 


This work 


^ 89 
O.oZ 


a 1 1 n 

41 1 .U 


4.0 


0O.0 


mC32 


This work 


0.10 


oo4.o 




7n 9 


0CI6-II 


This work 


5.95 


408.4 


3.15 


84.4 


W-carbon 


This work 


6.04 


391.8 


4.35 


83.1 




This work 


6.06 


392.6 


3.51 


82.7 


M- carbon 


PBE^P. 


5.97 


392.7 


3.60 




Bet C 4 


This work 


6.11 


393.4 


2.6 


82.0 


PBE22 


6.01 


411.4 


2.7 





a According to the method of A.O. Lyakhov and A.R. Oganov 33 




FIG. 3. Electronic band structures (top) and phonon disper- 
sion curves (bottom) for (from left to right) mC12, 0CI6-I, 
mC32 and oC 16-11 carbon phases. 



depleted in the experiments^ such that the experimental 
match is substantially the same for 0CI6-I and oC 16-11. 

Since superhard graphite is not synthetized from the 
gas phase, which would probably produce oC 16-11 as the 
only product due to its lowest enthalpy, in the real ex- 
periment much will depend on the nature of the starting 
graphitic material, and on the particular nucleation his- 
tory, which would favor one pattern in the phase growth 
phase. In this context, the overall stability of a particular 
structure is not the only parameter. The importance of 
this point of view has been recently pointed out 34 , and 
dedicated investigations are ongoing. 



Jj L l_J 



— OC16-I 
OC16-II 



LlLAj 



10 



■Hi ..J. 

12 14 



AjJ 

16 18 



29 



FIG. 5. Simulated XRD pattern for 0CI6-I and 0CI6-II car- 
bon phases (A = 0.3329A). The structural data are those of 
Table [D 
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In conclusion, we have presented four novel sp 3 carbon 
materials, derived from combining met adynamics and 
topology to achieve higher scan efficiency. Two struc- 
tures, 0CI6-I and oC 16-11 stand out for hardness and 
band gaps, and should be considered in assessing the na- 
ture of the product of graphite cold compression. 
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